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ABSTRACT 

This modeling effort seeks to improve the interlaminate bond strength of thermoplastic carbon 
composites produced by the in-situ automated tape placement (ATP) process. An existing 
high productivity model is extended to lower values of the Peclet number that correspond to 
the present operating conditions of the Langley ATP robot. (The Peclet number is the 
dimensionless ratio of inertial to diffusive heat transfer.) In sensitivity studies, all of the 
process and material parameters are individually varied. The model yields the corresponding 
variations in the effective bonding time (EBT) referred to the glass transition temperature. 
According to reptation theory, the interlaminate bond strength after wetting occurs is 
proportional to the one-fourth power of EBT. The model also computes the corresponding 
variations in the thermal input power (TIP) and the mass and volumetric process rates. 
Process studies show that a 10 percent increase in the consolidation length results in a 20 
percent increase in EBT and a 5 percent increase in TIP. A surprising result is that a 10 K 
decrease in the tooling temperature results in a 25 percent increase in EBT and an 8 percent 
increase in TIP. Material studies show that a 10 K decrease in glass transition temperature 
results in an 8 percent increase in EBT and a 8 percent decrease in TIP. A 20 K increase in 
polymer degradation temperature results in a 23 percent increase in EBT with no change in 
TIP. 
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1. INTRODUCTION 

Carbon fiber composites have greater strength and stiffness per unit mass than most materials 
used in aerospace applications. However, the fabrication of aircraft wings and fuselages from 
such composites is relatively expensive. In-situ automated tape placement (ATP) is a process 
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aimed at reducing this cost. Thermoplastic carbon composite tape — typically of thickness 
T = 1.27 x 10 4 m and width w = 0.0762 m — is laid down and thermally welded under 
pressure to a previous layer of ribbon to build up a laminated part. 

The in-situ ATP process has been modeled extensively, as reported in references 1-10 
However, a question that needs further attention is the following: What design changes in 
either the process or in the bulk properties of the composite would increase the interlaminate 
bond strength, decrease the thermal input power, or otherwise benefit the in-situ ATP 
process? To address this question, the authors will extend a thermal model given in reference 
10 and perform sensitivity studies. The Peclet number Pe (dimensionless ratio of inertial to 
diffusive heat transfer) in the model will be reduced by two orders of magnitude so that it 
corresponds to the present operating conditions of the Langley ATP robot. The model 
determines the effective bonding time t eff referred to the glass transition temperature 7’ , the 

thermal input power Q ) , and the mass and volumetric process rates M and V . According 
to reptation theory, as presented in reference 9, the interlaminar bond strength builds as l l[t ' 14 
after wetting occurs until full strength is reached. 

In the model, t eff , Q, , M , and V are functions of the both the process variables and the bulk 

properties of the composite. The process variables include the tape width w and thickness 
T , the laydown speed U , the consolidation length x h , and the temperature of the tooling at 
infinity 7A The bulk properties of the composite include its density p, specific heat at 
constant pressure C p , and longitudinal and transverse thermal conductivities K ] , and K 22 . 

These bulk properties can be modified by design to some extent. Hence, they are included in 
the sensitivity studies. 


2. MODEL 

The model uses an extended consolidation head of length x h , as shown in Figure 1. Heat is 
applied at the nip (x = 0, y = 0) . The thermal input power per unit depth is given by 
do = Qo ! w ■ This head configuration is similar to that reported in reference 1 1 ; however, all 
consolidation heads, including hard rollers, have finite consolidation lengths because of the 
pliability of heated composites. The consolidation head and the supporting tool are thermally 
conducting. Far away from the nip, the tool and head are maintained at temperature 7A The 
conductivities of the tool and head are both taken equal to the average thermal conductivity of 
the composite substrate. As shown in reference 10, this assumption, together with the 
assumption that Pe » 1 , enables the temperature field T(x,y ) to be determined by the 
classical two-dimensional solution for a line thermal source in an infinite, moving, anisotropic 
conductor with the boundary condition that T = at infinity. 

The thermal input power density q Q is constrained so that the temperature on the weld 
interface r(x,0) decreases passively to the glass transition temperature T, at the 
downstream end of the consolidation head, where x = x h . This constraint insures that the 
composite does not expand after consolidation pressure is released. 


The classical solution is then nondimensionalized by introducing the following dimensionless 
(primed) variables: 
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The classical solution for the dimensionless temperature T' is then given by (ref. 10) 


cxp(bx')K 0 fy(x' 2 +y! 2 ) l2 ~\ 

T (*' ,y' ) = 1 J 

V ' exp (b)K 0 (b) 

and that for the dimensionless thermal input power density q 0 ' by 


cxp(b)K tl (b) 
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where K 0 is the hyperbolic Bessel function of zero order that is bounded at infinity. The 
thermal input power (in W) is given by 


a = 


2mv(T t -T.yjK^Kj 12 
exp (b)K 0 (b) 


(3a) 


For b > 10 , this result can be approximated by 

a » 2w(T g - T^ P C p K 2 Ux h )' 2 (b > 10) (3b) 

3. PARAMETERIZATIONS 

The solution (2a) for dimensionless temperature T (x , y ) depends upon only one 
dimensionless parameter b . Isotherms in the substrate (y' < 0) are plotted in Fig. 2 for 
b = 10, 100, and 1000. The isotherms for y' > 0 are mirror reflections of those shown. 
Thermal input occurs at x' = y' = 0 . Pressure is applied by the consolidation head over the 
range (0 < x' < 1). According to this figure, the model’s use of passive cooling avoids the 
problem of temperature springback on the weld interface downstream of the consolidation 
head, as experienced by a model that uses active cooling (ref. 6). 



Isotherm depth of penetration into substrate. The isotherms shown in Fig. 2 can be 
characterized by their width AT (7" ,b) and depth of penetration into the substrate 

/ max (T' ,b) . The depth of penetration / max (7" ,b) is plotted in Fig. 3. High penetration is 

desirable because reheating previous layers under the consolidation head promotes bonding. 
Therefore, the increasingly rapid falloff of y' max (7" ,b) as T' increases for b < 1000 

represents a loss of bonding efficiency. In reference 10, b was restricted to equal or exceed 
1000. Here, new parameterizations allow b to be as small as 10. For b > 10, y' max (7' ,b) 

is approximated by the parameterization 


y (T',b) 

J max v 1 ' 


.6065 10 10 ^ 


f .6065 X 
b ln T ' 


( T' V 


H. 


m(b)\ 


lod - 




b m T 


(b > 10;1 < T < 10) (4a) 


where 


m{b) = 1.808(£> - 8.829) -027 

(4b) 

sib) = 1.576 1(£>+.4805)' 7821 

(4c) 

H a ( fi) = + tanh («/3)] 

(4d) 


Here H a {j5 ) is the hyperbolic tangent representation for the unit step function, as discussed 
in reference 12, where the index a controls the rise steepness near f = 0. The exact and 
parameterized curves for y' max versus T' are shown for b = 1 0 , 100, and 1000 in Fig. 4. 


Isotherm Width. Since the consolidation head applies pressure only on the interval 
(0 < .T< 1), any heating upstream of the nip does not contribute directly to bonding. Hence, 
the total isotherm width does not contribute to bonding. The part that does can be 
conservatively taken as the positive segment AT + ( 7" ,b) of the x’ - axis that is enclosed by 

the 7' -isotherm. This segment is plotted versus T' for various b in Fig 5. Again, the 
increasingly rapid falloff of AT + (T' . b) as T' increases for £><1000 indicates a loss of 
bonding efficiency. For b > 10 , these curves are approximated by the parameterization 
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The exact and parameterized curves for Ad + versus T are shown for b = 1 0 , 100, and 
1000 in Fig. 6. 


4. TEMPERATURE VERSUS TIME 

The accumulated time interval At(T) that the temperature at a material point on the weld 
interface is equal to or greater than T while pressure is being applied can be determined as 
follows: During the initial heating, the temperature at the point equals or exceeds T while 
under pressure for a period A x + (T) / U . As additional plies are added, the point can be 
reheated to temperature T or above. The total number of such heatings, including the initial 
heating, is given conservatively by y max (T) / z , where T is the ply thickness. Each time the 

point is reheated, the period when the temperature equals or exceeds T becomes shorter. 
This shortening can be accounted for by a factor n / 4 , which is the area ratio for an ellipse 
and its circumscribed rectangle. The resultant formula for A t(T) is given by 

At(T) = (6) 

4 U z 

Test computations with a wide range of isotherms and values for b confirm that (6) is a 
moderately conservative formula. It can be expressed in terms of the dimensionless variables 
by using (1), as follows 


K 

A «T') = j 


/ \l/2 , 

( x 2 

Uz 


v 22 


7 j-M,(ry m (r) 


Substitution of the parameterizations (4) and (5) then yields 
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where E(T' ,b) is the bonding efficiency discussed in section 3. It is given by 
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The bonding efficiency has the range (0 < E(T' ,b) < 1) and a contour plot of E(T' ,b) is 
shown in Fig. 7. 

For a particular composite and given process variables, equation (8) can be used to calculate 
Af(T')as T varies over the range (\ <T'< TJ ). Figure 8 shows a plot of T versus At 
for the baseline example of the Fangley ATP robot applying PEEK/carbon composite tape, as 



described in section 6. Note that T' is capped at 7/ to prevent degradation. Physically, 

this could be done by slightly defocusing the thermal input beam, as discussed in references 2 
and 5. Since the lower end of the interval At can always be taken as zero, Fig. 8 is also a 
plot of 7' versus t . 


5. EFFECTIVE BONDING TIME 

For a nonisothermal process, the effective bonding time t eff relative to an isothermal process 
at T g is given for a typical thermoplastic composite by (ref. 9 ) 




l 7 .44[7t/)-r ] 
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(9) 


where the integrand is the universal WLF expression and the integration includes all periods 
for which Tit) > T g . This equation is valid for temperatures Til) within about 100 K above 

T g . However, for PEEK, T d is 410 K above T g , and an Arrhenius type of temperature 

dependence applies. Such a temperature dependence is included in an extended WLF 
expression, as given by Ferry (ref. 13) 
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If written in terms of the dimensionless temperature T' , the extended WLF form for t eff 
becomes 
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When integrated numerically for T' versus t , as plotted in Fig. 8 for the baseline example 
described in section 6, equation (1 1) gives t eff = 7.33 x 10 13 s. Finally, the mass and 
volumetric process rates are given by 

M = pUxw (12) 


and 


V = Utw 


(13) 



6. SENSITIVITY STUDIES 


Baseline Example. Baseline process parameters, denoted by superscript 0, that 
approximate the present operating conditions of the Langley ATP robot are given by 

w° =0762 m 

T° = 1.27 x 10" 4 m 

U°=. 085 ms' 1 

x h ° = 1.27 x 10" 3 m 

T 0 = 293.16 K (14) 

Baseline material parameters for PEEK/carbon composite are given by (ref. 3 ) 

p° = 1560 kg m' 1 
C p ° = 1425 J kg 1 K 1 
Kj = 6 Wm 1 K 1 
K 22 ° = 72 W m 1 K 1 
T g ° = 413.16 K 

T d ° = 823.16 K (15) 

For these baseline process and material parameters, the model yields the baseline results 

b° = 20 

Q,° = 428.7 W 
t eff ° =7.328 xlO 13 s 
M° = 1.283 xlO" 3 kgs' 1 

V° = 8.227 x 10" 7 m 3 s' 1 (16) 


Sensitivity Studies. The computed results of individually varying the process parameters 
from their baseline values are summarized in Table 1. 



Table 1. Process Sensitivity Study 
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The computed results of individually varying the material parameters from their baseline 

values are summarized in Table 2. 






Table 2. Materials Sensitivity Study 
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7. DISCUSSION 



The sensitivity results for 

blb 0 , QJQ 0 , 

M I M° , and V I V° follow directly from (lc). 


(3a), (12), and (13). The results for t eff / t eff ° are more complex, since they follow from (11), 
where T as a function of t (or A t ) is given implicitly by (8), as illustrated for the baseline 
example by Fig. 8. Experience shows that the dominant contributions to t eff come from the 
7' values that are close to 7/ , where in the baseline example T d ' = 4.167 . 

Table 1. In Table 1, a 10 percent increase in x h causes a 20 percent increase in t eff , which 
is due to a 20 percent increase in A t(TJ ) . Fifteen percent of this increase comes from x h 3/2 
and 5 percent from an increase in efficiency E(T d ' ,b) that results from the 10 percent 
increase in b . The 10 percent decrease in U causes a 12 percent increase in t eff that is due 
to a 12 percent increase in At(T d ' ) . This increase results from a 17 percent increase due to 
U~ V2 and a 5 percent decrease in E(T d ,b) due to the 10 percent decrease in b. The 10 
percent decease in thickness T causes an 1 1 percent increase in t eff due to r 1 in (8a). The 
10 percent decrease in tape width w leaves t ejJ unchanged. A surprising result is that a 10 
K decrease in tooling temperature 7^ gives a 25 percent increase in t eff . The explanation is 
that T d is decreased by 6 percent, which gives a 25 percent increase in A t(T d ) from a 20 
percent increase in ( 7/ )“ 3 and a 5 percent increase in 7(7/,/?) . 



Table 2. In Table 2, a 10 percent decrease in either p or C p gives only a 1 percent 
increase in t eff because the 5 percent increase in A t(TJ ) due to p 1/2 or C~ 112 is 
counteracted by a 4 percent decrease in E(TJ ,b) due to the 10 percent decrease in b . The 
10 percent decrease in longitudinal conductivity K x , gives a 4 percent increase in t eff 
because the 1 1 percent increase in b gives a 4 percent increase in E(T d ' ,b), which results in 
a 4 percent increase in At(TJ ) . The 10 percent increase in K :1 gives a 5 percent increase in 
t eff simply because of the factor K 22 2 in A t(TJ ) . 

Up to this point, the sensitivity studies may be viewed as stretching the time axis of the 
integral in (11). In the subsequent studies, the amplitude of the integrand will also be varied. 

A 10 K decrease in T, gives an 8 percent increase in t eff , as follows: This decrease in 7' 

results in an 9 percent increase in TJ , which causes a 32 percent decrease in A t(TJ ) due to 
a 23 percent decrease in (7^')“ 3 and a 9 percent decrease in E(TJ,b). However, this 
decrease in At(TJ ) is more than offset by a large increase in the integrand at T = TJ to 
give t eff a net increase of 8 percent. Finally, a 20 K increase in T d results in a 23 percent 
increase in t eJf , as follows: This increase in T, causes a 4 percent increase in 7/ , which 
decreases A t{TJ ) by 14 percent because of a 10 percent decrease in ( 7/ )“ 3 and a 4 percent 
decrease in E(T d ',b). However, this decrease in A/ (7/ ) is more than offset by a large 
increase in the integrand at 7' = TJ to give t eff a net increase of 23 percent. 

The sensitivity studies show that t eff can be substantially increased without decreasing the 
productivity rates M and V by increasing the consolidation length x h , decreasing the 
tooling temperature 7^, or increasing the degradation temperature T d . Smaller increases in 
t eJJ can be obtained without decreasing M and V by decreasing the longitudinal 
conductivity K x , , increasing the transverse conductivity K 22 , or decreasing the glass 
transition temperature 7, • The thermal input power Q t can be reduced without decreasing 
M and V by decreasing the specific heat at constant pressure C or decreasing the glass 
transition temperature 7 . Decreasing the density p can also reduce Q } without reducing 
V , although M will be reduced. 

8. CONCLUSIONS 

An earlier thermal model for in- situ automated tape placement (ATP) has been extended to 
lower values of the Peclet number. This extension brings the model to the present operating 
Peclet numbers of the Langley ATP robot. Sensitivity studies were performed about these 
operating conditions with PEEK/carbon composite. The model computes the effective 
bonding time referred to the glass transition temperature, the thermal input power, and the 
mass and volumetric process rates as functions of the process and material variables. The 
process variables are the consolidation length, the laydown speed, the width and thickness of 
the tape, and the temperature of the tooling. The material variables are the density, the 
specific heat at constant pressure, the longitudinal and transverse thermal conductivities, and 
the glass transition and degradation temperatures. The computed thermal input power always 
causes the temperature on the weld interface to decrease passively to the glass transition 
temperature at the downstream end of the consolidation head. 



According to reptation theory, the interlaminar bond strength after wetting has occurred is 
proportional to the one-fourth power of the effective bonding time until full bonding is 
achieved. The sensitivity studies can be used as guides toward increasing the effective bonding 
time, decreasing the thermal input power, increasing the process rates, or otherwise optimizing 
the ATP process. A surprising result was that a 10 K decrease in the tooling temperature 
causes a 25 percent increase in the effective bonding time and an 8 percent increase in the 
thermal input power. The explanation is that reducing the tooling temperature also reduces 
the dimensionless degradation temperature, which increases the bonding efficiency. A 10 
percent increase in the consolidation length was found to give a 20 percent increase in the 
effective bonding time and a 5 percent increase in the thermal input power. A 20 K increase in 
the polymer degradation temperature gives a 23 percent increase in the effective bonding time 
with no change in the thermal input power. Other smaller, but significant changes in effective 
bonding time, thermal input power, and process rates were obtained as each of the process and 
material parameters were individually varied. 
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Figure 2. Concluded. 













Figure 7. Contour plot of efficiency F.(T‘ ,b ) , as given by (8b). 



Figures. Dimensionless temperature T' versus accumulated lime A/(D that the 
dimensionless temperature of a material point on the weld interface equals or exceeds T' , as 
given by (8), Abscissa Af can be taken as / for numerical evaluation of integral in ( 1 1). 






